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ABSTRACT 

Most theoretical investigations of particle acceleration during solar flares can- 
not be applied to observations for detailed study of the time evolution. We pro- 
pose a phenomenological model for turbulence evolution and stochastic particle 
acceleration that links observations to the energy release and particle acceleration 
through two coefficients characterising particle interactions with turbulent elec- 
tromagnetic fields. In the linear regime the particle distribution does not affect 
the turbulence energy cascade. It is shown that electron acceleration critically 
depends on the intensity of small-scale turbulence and an impulsive nonthermal 
component only appears near the peak of the gradually evolving turbulence in- 
tensity. The model naturally reproduces the soft-hard-soft pattern of hard X-ray 
pulses, and we attribute the observed change in fiux and spectral index correla- 
tion from the rise to decay phase of some pulses to changes in the background 
plasma. Detailed modelling of well-observed individual events will probe the 
energy release processes. 

Subject headings: Sun: fiares — Sun: particle emission — acceleration of particles 
— plasmas — turbulence 



1. INTRODUCTION 



Given the short particle acceleration timescales inferred from observations of solar fiares, 
it has been assumed that t he acceleration pr o cesses are decoupled from other processes 



in most theoretical models ( Miller et al. 1 119971 : lAschwanden I l2002l : iPetrosian fc Liu 1 12004 



Fletcher &: Hudson II2008I ). As a result, theoretical investigations of particle acceleration 



are aimed at meeting some basic observational requirements, such as timescales, energet- 
ics, and numbers. There are significant ambiguities in quantitative tests of specific models 
with observations due to free parameters introduced to characterise the related, not-well- 
understood processes. On the other hand, particle acceleration is just one important aspect 
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of the energy release processes, which are essentially multi-scale as evident from the rich 
time, space, and energy characteristics of flares. Direct modelling of the energy release, 
which can be related to magnetohydrodynamical simulations of large-scale processes, may 
better our understanding of flares. In the context of stochastic particle acceleration (SA) 
by turbulent electromagnetic fields, we introduce the concept of elementary energy release 
events characterised with an energy release rate and length scale and discuss its implications. 

The relation between the gradually varyin g thermal and bursty nonthermal emission 
component is a critical aspect of flare studies (IVeronig et al. II2005I ). Previous modelling 
focused on the correlation between the light-curve of nont hermal emission and derivative 
of the thermal emission light-curve (jPennis fc Zarro Ill993l ). This correlation has been at- 
tributed to chromospheric evaporation driven by a beam of nonthermal electrons that is 
injected at the top of flaring coronal loops, propagates to t he chromosph ere along magnetic 
field lines, and coUisionally heats the background plasma (IBrown Ill973l ). In this 'classical' 
electron beam model, most of the nonthermal ener gy is reprocessed to appear as therma l 
energy producing the dominant thermal signatures (INeupert Ill968l : lAntonucci et al. Ill984l ). 
However, this model does not address the origin of impulsive nonthermal electrons, and has 
difficulties in accounting impulsive soft X-ray emission observed from some footpoints and 
slower t han expected decay of t hermal emission after the nonthernial emission already dis- 
appears (IDennis fc Zarro Ill993l : iHudson et al. lll994l : lLi et al. 1119971 : 1 Jiang et al. II2006I ). By 
assuming a low-energy cutoff for a power-law or broken-power-law nonthermal electron dis- 
tribution, observations give a very poor constraint on the energetics since the beam power is 
dominated by low-energy electrons due to steep power-law slopes, and there are uncertainties 
in observationally constraining t his low cutoff energy due to the dominance of lo w-energy 
X-rays by a thermal component (IHolman et al. II2003I : ISaint-Hilaire fc Benz 1120051 ). Such a 
cutoff is also not expected in most theoretical mode ls hampering further investigations even 
for some well-observed flares (IVeronig et al. II2005I ). With the elementary energy release 
events, we propose a simple alternative to this electron beam model. The model is based 
on two facts that, 1) while low-energy particles couple strongly with each other through 
Coulomb collisions forming a thermal distribution, 2) high-energy particles must decouple 
from the th ermal background and obtairi their energy through interactions with electromag- 
netic fields ( iHamilton fc Petrosian Ill992l ). The key issues are then the dynamics of the fields 
and their interactions with charged particles. 

A significant amount of energy can be released during fiares, and the fiux of nonther- 
mal particles inferred from observations in combination with the source size measurement 
sometimes implies a no nthermal particle density possibly comparable to that of the coro- 
nal background plasma (iFletcher fc Hudson II2008I ). suggesting very efficient acceleration. A 
large fraction of the solar corona must be involved if the bulk of the energy and accelerated 
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particles are stored in the pre -flare coronal magnetic fields, as postulated in most theoretical 
models (iGrigis &: Benz II2006I ). On the other hand, the strong coronal magnetic field implies 
tiny gyro-radii for charged particles. An energy cascade from large to small scales appears to 
be inevitable. Turbulence is the most natural agent for energy cascade over a large dynam- 
i cal range, and the corresponding SA mod els have been widely used for solar fiare studies 
(iMiller et al~lll997l : IPetrosian fc Liull2004[ ). 

In these models a broad spectrum of particles is energised by interacting stochastically 
with a spectrum of electromagnetic fiuctuations over a broad range of spatial and temporal 
scales. One of the key features of the SA is that wave-particle interactions determine not 
only the energy gain rate of particles but also their spatial diffusion along magnetic field 
l ines. Consequently, the particle a cceleration is very sensitive to the turbulence intensity 
(IBenz II1977I : IPetrosian Sz Liu 1120041 ). which may explain the bursty behaviour of nonthermal 
emissions. Detailed studies of nonthermal hard X-rays (HXRs) also reveal spectral soft- 
hard-soft (SHS) evolu t ion of HXR pulses sugges ting independent electron acceleration events 
(IGrigis fc Benz 112004 : iBattaglia fc Benz 1120061 ). We propose a phenomenological model for 
the evolution of turbulence associated with elementary energy release events and use two 
coefficients to characterise electron interactions with the turbulence (Section [2]). The model 
can naturally fit HXR observations. With imaging spectroscopic observations of RHESSI, 
the model will allow us to extract the turbulence properties and its evolution for individual 
well-observed fiares, which constrain the wave-particle interactions and better our under- 
standing of the relation between the thermal and nonthermal emission components (Section 
[3]). Conclusions are drawn in Section |H Notations for different quantities are listed in Table 

m 



2. Elementary Energy Release Events and Stochastic Particle Acceleration 



(Benz 


I977I: 


Miller et al. 


1997 


Bvkov 


& Fleishman 


J2OO9) 


moc 



to an injection process, which is treated as an input independent of the turbulence evolution. 
Observations of solar flares, on the other hand, indicate that the large-scale energy release 
and particle acceleration are intimately connected. We aim at realising such a connection 
with a phenomenological SA model so that the model can be applied to individual events 
for detailed studies. 



We consider elementary energy release events, for which large scale eddies are assumed to 
be generated instantaneously upon the event trigger, which can correspond to a simple flare. 
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Table 1. Notations 



Quantities 


Symbol 


Typical Value 


Units 


Basic Quantities 


Magnetic Field 


B 


100-500 


Gauss 


Density 


n 


lOMQio 


cm"^ 


Flaring Region Length 


h 


lOMQio 


cm 


Loop Cross Section 


A 


lO^MO^^ 


cm^ 


Temperature 


T 


lOMO^ 


K 


Energy Release Scale 


le 


lOMo^ 


cm 


Energy Release Level 


bo 


< 1 


1 


Derived Quantities 


Alfven Speed 




108-109 


cm s ^ 


Turbulence Intensity 


b 


<bo 


1 


Eddy Speed 


Ve 


b 


Va 


Transition Energy 


Et 


1 




Particle Speed 


V 


> 109 


cm s--*^ 


Wave Transit Time 


To = le/vA 


0.01-10 


s 


Eddy Turnover Time 


Te = le/Ve 




To 


Rise Time 


Tr = To/bo 


bo' 


To 


Decay Time 


Td = TeVA/Ve 




To 


Particle Transit Time 


Tt{v) = le/v 


< 0.01-1 


S 


Scattering Time 


Tsciv) 


Ssb-' 


Tt 


Acceleration Time 


rac{v) 


Sab-' 


Tt 


Escape Time 


resc{v) 


72 -2 -2 
'O " ^sc 


Tsc 
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or one HXR pulse within a large complex flare. These eddies are characterised by a generation 
scale le, and speed Veo = &o^a- The particle acceleration from the background plasma is 
determined by small-scale turbulence characterised by the amplitude of the corresponding 
magnetic field fluctuations bB. The growth of b is driven by large-scale eddies with growth 
rate = Veo/h = bo/ro, where tq = h/vA is the transit time of Alfven waves through Ig. 
The turbulence starts to decay once b reaches 60, and we adopt the Kraichnan phenomenology 
with the decay time given by = IpVa/v"^ = TpVa/ vp, where the eddy speed = bvA and 
the eddy turnover time Xg = /g/fe (IKraichnan Ill965l ). Then we have 



b/b 



bo/ro 
= -bVro 



for the turbulence rise phase , 
for the turbulence decay phase . 



This is the pr oposed basic equati on for the turbulence evolution. Considering the variety of 
flare triggers (lAschwanden II2002I ). the turbulence evolution can be much more complicated. 
These equations can be modified wherever there are sufficient observational or theoretical 
justifications. 




Fig. 1. — Evolution of elementary energy release events. The time unit is tq/A and tq/A^ 
for the rise and decay phase, respectively. Upper: the thick and thin lines are for b^ = A and 
A/2^/^, respectively. The latter lines are shifted to the right by ro/2A^ so that the decay 
phase overlaps with the former. is for Eq = lOksT and normalised to the peak value for 
60 = A. Lower: the dotted line is for A^E. The solid line and normalisations of A^ and E are 
the same as the upper panel. The dashed lines are the same as the thick-dashed line in the 
upper panel except that from high to low the corresponding electron temperatures are 0.1, 
5, and 10 times lower, respectively. 



We consider the solution, where b = 
uu-^ _ u / exp(t/r,) 



60 at t = 0: 

for t < the turbulence rise phase , 
for t > the turbulence decay phase . 
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The solid lines in Figure [T] show the time evolution of 6/A for bo = A (thick) and 
(thin), where A < 1 is an upper limit for bo, below which the turbulence evolution is not 
significantly affected by nonthermal particles (see discussions near the end of this section). 
The corresponding total amount of energy per unit volume dissipated through small scale 
turbulence is given by 

m) - 2ov^ t b\-'dt' -avH'i exp(4t/r,) for t < , 

L(t) - 2pv^ J 6 r, dt - pv^b, I ^^^2 + 1 - (1 + 2t6g/ro)-i for t > , 

where p is the mass density. It is indicated by the dot-dashed lines in the Figure. The 
total energy dissipated through turbulence is (bo/2 + l)pv\bl with the first and second term 
resulted from the rise and decay phase, respectively. Less than 1/3 of the total released 
energy is dissipated in the turbulence rise phase. Most of the energy is dissipated in the 
much longer decay phase. The rise of E(t) approximately represents the rise of thermal 
emission. 

We characterise the decoupling and acceleration of individual particl es from the thermal 



background by small scale turbulence with an acceleration timescale (IMiller et al. I Il997 



Petrosian fc Liu 1120041 ) 



rac = SaTtb-\ (4) 

where Sa is a dimensionless coefficient describing the acceleration by waves. The energy 
loss time ti{E) of high-energy particles through Coulomb collisions with a low temperature 
background is proportional to E^^'^/n, where E ^ ksT is the particle kinetic energy, n and 
T are the background particle density and temperature, and k-Q is the Boltzmann constant. 
The particle distribution is nonthermal at energies, where Tac{E) < ti{E). The transition 
from the non-relativistic thermal to nonthermal component occurs at 

Et = 2[7rlnAe^nr,,/mi/Y^^ , (5) 
where e, m are the particle charge and mass, respectively , and InA ~ 20 for plasmas in the 



solar corona, i.e., Tac{Et) = Ti(Et) fjPetrosian &: Liu II2004I ). 



The most natural way to produce a power-law high-energy particle distribution with 
an adjustable spectral index is in terms of a particle loss process as proposed originally by 



Fermi I (119491 ) for cosmic rays. For solar flares, this particle loss is due to escape from the 



energy release site through spatial diffusion. We use the scattering time 

Tsc = SsTtb'"^ (6) 

to characterise this spatial diffusion, where Sg is the second coefficient describing the particle 
scattering by waves. The corresponding particle escape time 
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where Iq > le is the length of the energy release site. The electron acceleration timescale is 
much shorter than the flare duration. One may consider the steady-state solution. To have 
a power-law distribution as commonly adopted to model the observed nonthermal emissions, 
''"ac and Tgsc need to have the same energy dependence, i.e ., TacTsc l/v'^. and Tac may 



have different energy dependence (jPetrosian fc Liu II2004I ). For simplicity, we assume that 



their energy dependence is the same. For non-relativistic particles, as considered in the 
paper, v cx i?^/^, therefore Tac oc and both 5"^ and Ss are independent of v. 

In the high-energy range, where Tac ^ Ti, the kinetic equation for nonthermal particle 
distribution f{p) is given by 

— = ——p^r-^—f - — + Q (8) 

dt p2 dp ^"^ dp Tesc ' 

where p is the particle momentum, and Q, a source term, exists at low energies. In the 
steady- state. 

The corresponding energy distribution N{E) oc iv;^i/2-(i+f"ac/4Tesc)^''^ ^ g^^^j ^i^q of escaping 
particles F{E) oc N{E)/Tesc oc E-(i+^-/4re.,)V2 _ 

convergence of energy flux carried 
away by nonthermal particles, the index of F{E) must be greater than 2, i.e., 

Tac > 12W , < SaSsile/lo)Vl2 = . (10) 

For b > A, one must consider the relativ istic effect and/or the effect of nonthermal particles 



on the turbulence cascade and damping (IBykov fc Fleishman II2009I ). We consider the linear 



regime, where nonthermal particles do not affect the turbulence cascade. A is an upper limit 
for the linear model to be valid. 



3. Application to Solar Flares 

Electron acceleration during flares is better observationally constrained than that for 
ions. We next apply the model to electron acceleration in flaring loops with ions treated 
as a background of positive charges, which only changes through large scale hydrodynamic 
processes. □ Then the length of the energy release site Iq has a lower limit of the size of coronal 



^ Other assumptions of the energy dependence of Tac a-nd r^c wiU lead to different energy dependence of r, 



which determines the flux of escaping particles and can be constrained by observations ([Battaglia fc Benz 



20061 ) . For example, with the standard quasi-linear theory, Tacj'f'sc oc (w/wa)^- It can be shown that t^sc is 



independent of E and 5a oc 1/55 oc u. 

^This is an assumption. Ion acceleration can be considered in appropriate theoretical and/or observational 
contexts. 
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loopt op sources (IXu et al. II2008I ) and an upper limit of the length of flaring loops (iLiu et al 
20061 ). Although high-energy electrons must escape from the acceleration site to produce a 
power-law distribution with the above SA model, charge neutrality of the acceleration site 
requires the total number of electrons remains the same as positive charges at Nq = IqAu 
(e.g., due to a return current), which gives the normalisation for the electron distribution 
at the acceleration site N{E). One of the key purposes of the paper is to introduce a 
simple phenomenological model to extract the energy release rate through turbulence from 
observations. Instead of solving the full kinetic equation for the electron distribution over 
the whole ene rgy, i.e., from the Cou l omb collision dominat ed low-energies to coUisionless 
high-energies JPetrosian fc Liulbooi ICallowav et al. Iboosh . we assume that the electron 



distribution is thermal below Et, a power law above it, and continuous ai E = Et. Then 
N{E) = NqqIE), where g{E) is the normalised distribution function, i.e., f g{E)dE = 1 and 



g^{Et/k^T,5) f ^1/2 eM-E/k^T) for E < Et , 

(k^T)^/^ \ Ey'exp{-Et/kBT){E/Etr' for E > Et , 



where 

5 = 1/2 + (1 + r,,/4w)^/2 ^ (^2) 
gQ(r, S) = [r=^/' exp(-r)/(5 - 1) - r^/^ exp(-r) + 7r^/hTf{r^/^)/2]-^ , (13) 

and erf is the error function. Et = {2n\nASJeny/^e^ /b , 5 = 1/2 + (1 + SA^/b'^y^^ ,r = 
Et/k^T = R{A/b) with R = {2T:\nKSalenf''^e^/k^TA. 

N is described with four parameters: Nq, T, Et, and S. For a given energy release 
event discussed in Section |2] with bo given in units of A, the evolution of 6 is determined 
as indicated by the dotted lines in Figure [H For bo = A, 6 reaches 5/2 at the peak of b. 
Nq and T can be measured observationally and vary gradually as S on a timescale much 
longer than Tac- We assume that they do not change during these events. For Et = k^T at 
b = A (the fiducial model), the evolution of N at Eq = lOkBT (normalised to the peak value) 
is indicated by the dashed line. Because 6 has a strong dependence on b, the nonthermal 
electron density at high energies is very sensitive to b. A prominent peak of nonthermal 
electron numbers only appears near the peak of b. The thin lines show another event with 
about 2 times lower released energy. As expected, the electron distribution becomes softer 
than the previous model with a minimum of 5 at 1/2 + 13^/^ ~ 4.1. The peak density of 
nonthermal electrons is smaller by more than one order of magnitude. The decrease of the 
peak density is more prominent at even higher energies due to the softer spectrum. Therefore 
strong and hard nonthermal emissions are expected only for strong events with fog close to 
A, and even with a relatively gradual evolution of the turbulence intensity, an impulsive 
nonthermal component appears at the peak of the energy dissipation. This may explain 
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the bursty nature of the nonthermal component as compared with the thermal component. 
The dependence of N{Eq) on T is shown in the lower panel. There are fewer nonthermal 
electrons for lower T. The evolution of T and A^^o will affect the quantitative details but not 
the impulsive nature of the nonthermal density, which is mostly determined by S{t). 




-12 -10 -8 -6 -4 -2 2 

Ln(N), Ln(NS} 



Fig. 2. — Correlation between nonthermal electron spectral index and density. The solid 
line is for the fiducial model. The other lines correspond to those in the lower panel of Figure 

m 



One of the most importan t observations o f nontherma l emission is the spectral SHS 
evolution of HXR pulses, which iGrigis fc Benz I (12004 . |2006| ) attributed to elementary elec- 
tron acceleration events. The spectral index and nonthermal electron density evolution of 
our model naturally reproduces such a pattern. The thick solid line in Figure [2] shows the 
correlation between S and N{Eq) of the fiducial model. The correlation at a given energy is 
the same for the rise and decay phase in agreement with observation s of some HXR pulses . 



For most pulses, this correlation changes from the rise to decay phase (IGrigis fc Benz 



2004). 



This can be caused by chromospheric evaporation and/or plasma heating due to the energy 
release S, which changes the density and temperature of the background plasma. The dashed 
lines show how the correlation changes with T and the dotted line shows the correlation be- 
tween S and NH. The chromospheric evaporation processes must be modelled properly to 
study this correlation. 



4. Conclusions 



Yohkoh and RHESSI observations of solar flares have revealed several phenomena chal- 
lenging the classical electron beam model: impulsive soft X-ray emission from footpoint 
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sources ( iHudson et al. Ill994l : iHannah et a. 



after the impulsive HXR phase (ILi et al 



1997 



20081) . continuous heat ing of the coronal source 



Jiang et al. 1 120061), and more recently in- 



jection of electrons over a large coronal looptop region (jXu et al. I l2008l ). All these point 
to the SA model, where electrons are accelerated by magnetised turbulence in flaring loops 
( jPetrosian fc Liu 112004 ). The acceleration site can be as compact as the observed coronal 
looptop sources and may also extend over the whole loops if the turbulence is transp orted 
along the loops quickly, e.g., in the form of plasma waves (IFletcher fc Hudson II2008I ). The 
latter may explain impulsive soft X-ray emission from footpoints as electrons are injected at 
the footpoints directly. In this paper, we propose a simple phenomenological SA model and 
show that the nonthermal electron density in the acceleration site is very sensitive to the 
turbulence intensity. It therefore provides a mechanism to produce impulsive emissions even 
with a relatively gradual energy release process, which may account for the continuous heat- 
ing inferred from thermal source evolution after the impulsive phase. With a simple model 
for the turbulence evolution, it also reproduces spectral SHS evolution of HXR pulses. De- 
tailed modelling of the plasma heating and chromospheric evaporation is needed to quantify 
the flux and spectral index correlation. 

Solar flares are multi-scale phenomena in terms of not only the size and duration but 
also the amount of energy released by each flare. The physical processes involved in the large- 
scale energy release process can be scale dependent, which in combination with the variety 
of initial an d boundary conditio ns in the solar atmosphere is expected to lead to very rich 
appearance (lAschwanden II2002I ). However, the microscopic scale processes of plasma heat- 
ing and particle acceleration should be mostly determined by properties of the background 
plasma in terms of temperature, density, and large scale magnetic field. Besides Coulomb 
collisions among these particles, particle interactions with the electromagnetic field fluctu- 
ations during the energy release process may be parameterized with an acceleration and 
scattering coefficient. We show in the paper that these two coefficients, in combination with 
the turbulence intensity and generation scale, should be used to determine the electron dis- 
tribution in flaring loops and applied to individual flares. The model may not only overcome 
the observational challenges to the classical electron beam model but also address the in- 
triguing problem of the low-energy cutoff. The two coefficients should not vary significantly. 
The turbulence intensity and release scale le are determined by the macroscopic processes of 
flare triggers. Detailed analyses of flare X-ray emissions will lead to quantitative constraints 
on the wave-particle interactions and turbulence evolution and may also help to understand 
the related large-scale processes. 



This work is supported by the EU's SOLAIRE Research and Training Network at the 
University of Glasgow (MTRN-CT-2006-035484) and by Rolling Grant ST/F002637/1 from 



-li- 



the UK's Science and Technology Facihties Council. We thank Hugh Hudson, Eduard Kon- 
tar, John Brown, and Alec MacKinnon for helpful comments. 

REFERENCES 

Antonucci, E., Gabriel, A. H., & Dennis, B. R. 1984, ApJ, 287, 917 

Aschwanden, M. J. 2002, Space Science Reviews, 101, 1 

Battagha, M., & Benz, A. O. 2006, A&A, 456, 751 

Benz, A. O. 1977, ApJ, 211, 270 

Brown, J. C. 1973, Solar Phys. 31, 143 

Bykov, A. M., & Fleishman, G. D. 2009, ApJ, 692, L45 

Dennis, B. R., & Zarro, D. M. 1993, Solar. Phys., 146, 177 

Fermi, E. 1949, Phys. Rev. 75, 1169 

Fletcher, L., & Hudson, H. S. 2008, ApJ, 675, 1645 

Galloway, R. K., MacKinnon, A. L., Kontar, E. P., & Helander, P. 2005, A&A, 438, 1107 
Grigis, P. C., & Benz, A. O. 2004, A&A, 426, 1093 
Grigis, P. C., & Benz, A. O. 2006, A&A, 458, 641 
Hamilton, R. J., & Petrosian, V. 1992, ApJ, 398, 350 

Hannah, I. G., Krucker, S., Hudson, H. S., Christe, S., & Lin, P. R. 2008, A&A, 481, L45 
Holman, G. D., Sui, L. H., Schwartz, R. A., & Emslie, A. G. 2003, ApJ, 595, L97 
Hudson, H., S., et al. 1994, ApJ, 422, L25 

Jiang, Y. W., Liu, S., Liu, W., & Petrosian, V. 2006, ApJ, 638, 1140 

Kraichnan, R. H. 1965, Phys. of Fluids, 8, 1385 

Li, P., McTiernan, J. M., & Emslie, A. G. 1997, ApJ, 491, 395 

Liu, W., Liu, S., Jiang, Y. W., & Petrosian, V. 2006, ApJ, 649, 1124 



- 12 - 



Miller, J. A., et al. 1997, J. Geophys. Res. 102, 14631 

Neupert, W. M. 1968, ApJ, 153, L59 

Petrosian, V., & Liu, S. 2004, ApJ, 10, 550 

Saint-Hilaire, P., & Bcnz, A. O. 2005, A&A, 435, 743 

Veronig, A. M., et al. 2005, ApJ, 621, 482 

Xu, Y., Emslie, A. G., & Hurford, G. J. 2008, ApJ, 673, 576 



This preprint was prepared with the A AS IM^jX macros v5.2. 



